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TWO-DIMENSIONAL UNSTEADY LIFT PROBLEMS IN SUPERSONIC FLIGHT 


By Max. A. Heaslet and Harvard Lomax 


SUMMARY 

The variation of pressure distribution is calculated for a two- 
dimensional supersonic airfoil either experiencing a sudden 
angle-of -attack change or entering a sharp-edged gust. From 
these pressure distributions the indicial lift functions applicable 
to unsteady lift problems are determined for the two cases. 

Results are presented which permit the determination of 
maximum, increment in lift coefficient attained by an unrestrained 
airfoil during its flight through a gust. As an application of 
these results, the minimum altitude for safe flight through a 
specific gust is calculated for a particular supersonic wing of 
given strength and wing loading. 

INTRODUCTION 

The stud} 7 of the unsteady lift of wings in an incompressible 
medium has been developed along two different lines. In 
reference 1, R. T. Jones introduced the concejit of indicial lift 
functions for wings of finite aspect ratio and, using as a basis 
the work of Wagner (reference 2) on the two-dimensional 
potential theory of airfoils in nonuniform motion, has shown 
how the calculation of lift under various conditions of motion 
can be effected. In reference 3, Tlieodorsen considered non- 
steady motion in its relation to the general theory of aero- 
dynamic instability and the determination of the aerody- 
namic forces on harmonically oscillating airfoils. This latter 
approach has been extended to include high-speed problems 
and in two recent papers Garrick and Rubinow (references 
4 and 5) have given results on flutter and oscillating air-force 
calculations for wings in supersonic flow. 

The present report employs the method of attack intro- 
duced by Jones and considers the case of a two-dimensional 
airfoil moving supersonically in an arbitrary manner, pro- 
vided the assumptions of small perturbation theory are satis- 
fied. The principal contribution lies in the determination of 
indicial pressure distributions which are readily calculated in 
supersonic motion and from which indicial lift, drag, and 
pitching moments may be computed. From these results 
the indicial lift functions are calculated explicitly. The 
methods used to find the pressure distributions also afford 
considerable insight into the same problem for airfoils at 
subsonic speeds. 

As an application of the analysis, the results are applied to 
the special case of an unrestrained airfoil entering a sharp- 
edged gust. The resultant forces are found to be comparable 
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in magnitude, for Mach numbers in the neighborhood of 1.3, 
to those given in reference 1 for subsonic incompressible flow. 

SYMBOLS 

A acceleration factor 

a speed of sound 

c chord length 

C L wing lift coefficient 

c, section lift coefficient 

Ci a section indicial lift coefficient for angle-of-attack 
change 

c, s section indicial lift coefficient for wing entering gust 
F net lifting force on wing 

g acceleration of gravity 

M fi'ee-stream Mach number 

m mass 

A p difference in pressures between lower and upper 

surfaces of airfoil 

q frec-stream dynamic pressure ^ Po Fo 2 J 

s distance measured in half-chord lengths 

S area of wing 

p perturbation static pressure 

t' , t' time in seconds 

t transformed time variable (See equations (7).) 

T t/c 

u perturbation velocity component in x direction 
V„ free-stream velocity 

w perturbation velocity component in 2 direction 
w 0 z component of velocity of gust 

W weight of wing 

x,z Cartesian coordinates 

a angle of attack 

p wing-density parameter (2 m/pMSc) 

p perturbation density 

p a free-stream density 

<r area over which surface integral is evaluated 

I> perturbation velocity potential 

SUBSCRIPTS 

o free-stream conditions 

! variable of integration 

u upper surface 

f lower surface 
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ANALYSIS 

DERIVATION OF BASIC DIFFERENTIAL EQUATION 

The pressure distribution over an airfoil in a compressible 
medium is obtainable from the solution of a boundary-value 
problem associated with a particular second-order partial 
differential equation. The derivation of the linearized form 
of this equation, obtained under the assumptions of small 
perturbation theory, proceeds as follows : 

Let u,w be perturbation velocity components parallel, 
respectively, to the Cartesian axes x', z' and denote by p 
perturbation pressure, by p perturbation density, and by a 
the velocity of sound. Then if t' denotes time, V a is the 
constant free-stream velocity, and p 0 is the constant free- 
stream density, the linearized Eulerian equations are 

bu . v bu &P ' 

bt' + °dx' p 0 bx' 

r (i) 

bw T/ dw_ 1 bp 
bt °bx'~ pods' „ 

The linearized equations of continuity and state are, respec- 
tively, 



After the introduction of the perturbation velocity potential 
<f> integration of equations (1) yields 

§P + Fo ^ =_ p L „ ?, + COnsfcanfc (4) 

while equations (2) and (3) give the expression 




V 

Po a o 2 



( 5 ) 


The combination of equations (4) and (5) leads to the desired 
partial differential equation 


( 1 - 




d 2 $ 

d? 


2 M d 2 $ 1 d 2 4> 

a B bx'bt' a 0 2 bt' 2 ~ 


(6) 


where M is the free-stream Mach number. 

Equation (6) can be reduced to the normalized form of the 
two-dimenstional wave equation of mathematical physics by 
means of the transformation 

x—x'—Ma 0 t' \ 

z=z' > (7) 

t=a 0 t' ) 


In these variables the equation is written 

d 2 4> d 2 $ d 2 <f> 

W bx 2 dF -0 


( 8 ) 


In accordance with the assumptions underlying the deriva- 
tion of equation (8), its application to problems in airfoil 
theory is, of course, limited to cases where the induced ve- 
locities are small compared to the free-stream velocity and the 
effects of viscosity do not alter the results of the potential 
flow solution. 

The rectangular coordinate system associated with equa- 
tion (6) is fixed in the wing which is, in turn, immersed in a 
free stream of velocity V 0 directed along the positive x axis. 
The transformations introduced in equation (7) fixes the x, z 
coordinate system in space so that the airfoil moves in the 
negative x direction and the free-stream velocity is zero. A 
distortion of the time axis is also involved so that the differ- 
ential equation appears in canonical form. 

SOLUTION FOR GIVEN BOUNDARY CONDITIONS 

The boundary conditions which are to be satisfied have the 
same property as those encountered in steady-state thin- 
airfoil theory; that is, the prescribed data are given in the 
z=0 plane. The particular problems with which this report 
deals are those of finding pressure distributions over a flat 
plate. Thus, w will be specified over a portion of the 2=0 
plane and, elsewhere throughout the plane, loading must 
be zero. 

Two boundary-value problems are to be considered: first, 
the case of an airfoil either starting from rest at a given angle 
of attack a or experiencing along the entire chord a change a 
in stream direction without a pitching motion, and second, 
the case of a constrained wing entering a sharp-edged gust 
with a vertical velocity w 0 . In the former case the motion is 
that of an airfoil suddenly sinking without rotation. These 
boundary conditions are more readily pictured with the aid 
of figure 1. Figure 1 (a) shows the conditions which must 
be satisfied in order to solve the angle-of-attack problem. 
The trace of the leading edge of the wing traverses the line 
x= — Mt, while the trailing edge lies on x=c—Mt where c is 
chord length. The region bounded by these lines and the 
line t = 0 is the region occupied by the airfoil as time passes. 
Since the axes are fixed and the airfoil moves in the negative 
x direction, the velocity at which the airfoil travels deter- 
mines the inclination of the leading- and trailing-edge traces. 
Over the “area” occupied by the airfoil in the x, t plane, w 
must equal— V 0 a and elsewhere no jump in pressure can occur. 
The gust problem (fig. 1 (b)) does not differ essentially from 
the previous problem except that here the region over which 
the modification of w is effective is not entirely the region 
occupied by the airfoil but rather the region occupied simul- 
taneously by the airfoil and the gust. Fixing, for conven- 
ience, the edge of the gust along the t axis, this axis will form 
the right-hand boundary of the region over which w= — w 0 . 

Since the partial differential equation is linear and the 
solutions are therefore additive, these boundary conditions 
clearly should fit the following physical event: A wing of 
trapezoidal plan form, indicated in figure 2, flies at a steady 
lift and angle of attack prior to <=0; at t— 0 the wing either 
experiences a change in angle of attack a with no pitching 
motion or enters a sharp-edged gust of constant vertical 
velocity w 0 , the gust extending from its edge to all negative 
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(a) Supersonic wing after a sudden angle-of-attack change. 

(b) Supersonic wing entering a sharp-edged gust atx=0. 

Figure 1. — Boundary conditions for supersonic wings. 


values of x. The wing in each case is then restrained so 
that, relative to the original wind vector, the wing remains 
at an angle of attack a or, in the gust case, continues fixed 
at the same angle of attack. 

The solution to similar boundary-value problems has been 
discussed at length in reference 6. In that report the 
development was adapted to the case where equation (8) 
represented the steady-state equation for wing problems in 
three dimensions and the characteristic cones of the equation 
had the immediate physical interpretation of Mach cones or 
infinitesimal shock disturbances. Green’s theorem was 
applied to solve the boundary-value problems involved and 
it was shown that the solution obtained could be interpreted 



as a surface distribution of sources and doublets. Finally, 
the difficulties arising in the discussion of the singularities 
on the characteristic cones and the integration of the super- 
sonic doublets were overcome through the introduction of an 
integration technique which involved using the “finite part” 
of the given integrals. 

The methods and conclusions obtained in reference 6 can 
be adapted immediately to the problems discussed herein for 
the mathematical reasoning remains almost identical. The 
physical interpretations of the two cases must, of course, be 
modified. Thus, the characteristic cones of the differential 
equations, traces of which are shown in figures 1 and 3, are 
no longer the well-known Mach cones; rather, they represent 
the distance to which a disturbance occurring at a point 
fixed by the apex of the cone will travel in the time t. 
Despite the fact that such physical interpretations are 
undeniably useful in understanding and applying the results, 
the solution of the basic differential equation for the bound- 
ary values involved is quite independent of these material 



Figure 3.— Sketch indicating locations of regions A, B, and C and time intervals used 

in analysis. 


dissimilarities. It follows that the methods developed and 
discussed in detail in reference 6 can be applied directly to 
the given problems with only minor changes in notation. 

With reference to figure 1, it is evident that for supersonic 
flight the air ahead of the wing is unaffected by the approach 
of the wing and, further, that the induced velocities on the 
upper surface of the airfoil are independent of the shape of 
the lower surface. Consequently, the pressure distribution 
will be found on the upper surface, as if the airfoil section 
were symmetrical, and then, for the flat plate, the pressure 
distribution on the lower surface will be equal in magnitude 
and opposite in sign. As in reference 6, the solution to such 
a problem can be obtained from a distribution of sources. 

In the actual computation of the pressures over a section 
traveling at supersonic speeds, certain regions are con- 
veniently defined. These regions depend on the relative 
slope of the traces of the leading and trailing edges and the 
trace of the characteristic cone in the r, t plane (fig. 3) . The 
perturbation velocity potential is given by the formula 



dx-jdti 

—t 1 ) 2 —(z—x 1 ) 2 —z 2 


( 9 ) 
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where <r is the area in the x, t plane of the region occupied by 
the wing section and bounded by the trace of the forecone 
from the point x, z, t. If p u and p t denote, respectively, 
pressures on the upper and lower surface of the wing 


A p_p t —pu 4 v d<b\_ 4 

2 ~ 2 ~ V 0 2 W + °bx')~ V 0 M bt 


( 10 ) 


Using equations (9) and (10), direct calculation shows that 
for a sudden angle-of-attack change, the following relations 
hold 

Region A (between lines x = — Mt, x=—t, and x=c — Mt) 


A p 4a 

_ <Z ~VM 2 - 1 


(11a) 


It follows that in this region the steady-state Ackeret-type 
load distribution has been attained. 

Region B (between lines x=—t, x—t, and x—c—Mt ) 


Aj o 4a 

q = \ M- I 




1 Mx-f-t 

- arc cos — , ,, 
7T X + Mt 


VM 2 - 1 

ttM 


(f -fare sin f)] 


(lib) 


Region C (between lines x=t, t= 0, and x=c—Mt) 


Ap_4a 

q~M 


(11c) 


The result obtained for Region C is of particular interest 
since it holds for airfoils at subsonic as well as supersonic 
speeds. Moreover, the mechanics of the interaction be- 
tween the airfoil and the fluid are such that other methods of 
derivation, furnishing added insight into the nature of the 
phenomenon, may be developed. Consider a flat plate of 
infinite aspect ratio flying at a velocity V 0 either greater than 
or less than the velocity of sound a a in the undisturbed air. 
The airfoil is assumed to undergo a change in its motion at 
the time t — 0 so that subsequent to this time it has in- 
creased its angle of attack by the amount a. It follows that 
the sudden increment in lift can be calculated from a knowl- 
edge of the induced effects on the air produced by an added 
vertical velocity of the plate equal to — V 0 a. As a result of 
this vertical motion two plane Rayleigh waves will emanate 
from the plate, a compression wave from the lower surface 
and an expansion wave from the upper surface. The 
velocities of the wave fronts are equal to a 0 while the induced 
velocities in the waves are equal in magnitude to U 0 a. 

The lift on the wing can be determined in two ways: from 
impulse relations and from energy considerations. In the 
former case, assume that the forces per unit span on the 
upper and lower surfaces are j u and respectively. After 
an elapse of time At the wave fronts have advanced a distance 
a 0 At and each includes a mass of air equal to p 0 ca 0 At per unit 
span. From Newton’s second law of motion 

(fi—fu)At=(p 0 ca 0 At) (2 V 0 a) 


and, converting to lift coefficient, 


ji—fu 

2 Pof'Ye 


40,00: 4a 

~W~M 


Since the force is distributed uniformly along the chord this 
result is the equivalent of equation (11c). 

In the development of the theory of plane waves of small 
amplitude (see, e. g., reference 7) Rayleigh and Lamb have 
shown that the energy in a wave is divided equally into 
kinetic and potential energy. Denote kinetic energy by 


1 

2 p0 


| jf* J w 2 dxdydz 


where w is the perturbation velocity within the wave. Since 
the energy induced in the wave must result from work done 
on the plate, it can be seen that 

(W«)V 0 aA* = ! po j o V n Wcdx 
or 

(Ji—fu) At=2p 0 V 0 aca 0 At 


This equation is in agreement with the one obtained pre- 
viously. 

As a consequence of equation (lie) it follows that the 
starting lift coefficient of an airfoil is equal to 4a /M for both 
subsonic and supersonic flight. The magnitude of c t thus 
increases as M becomes smaller and for incompressible 
theory, where the velocity of sound is indefinitely large, 
must necessarily become infinite. This fact was known 
previously along with the understanding that the indicial 
lift function experiences an infinite discontinuity at t= 0. 
(See Wagner’s curve in a subsequent figure.) For values of 
M other than zero the starting lift is finite and a continuous 
lift function results. 

For a sudden gust with vertical velocity w 0 the following 
expressions can be found for the corresponding regions 

Region A 


Region B 

A p 4w 0 

2 U v 'M 2 -7 


A p 

4w o 

Mx-\-t 


“ XT '-*-*<* ^ 

IT Vo -x M 2 — 1 

x+Mt 


(12a) 

(12b) 


Region C 



(12c) 


APPLICATIONS AND DISCUSSIONS 


DISCUSSION OF LOAD DISTRIBUTIONS 


Figure 4 shows the variation of the loading on a section 
which, while traveling at supersonic speed, is suddenly de- 
flected to a new angle of attack. The loading varies accord- 
ing to equation (11) for each of the three regions A, B, and C 
of figure 3. At t— 0 the pressure is discontinuous, jumping 
from its original value, just before the sudden deflection in 
angle of attack, to 4 a/M just after the deflection. Figure 3 
shows, however, that the initial load distribution is modified 
over the forward portion of the section as time increases 

c 

since regions A and B must be considered. For 
the wing lies entirely in region A so that 


A p _ 4 a 

2 _ tJM 2 — 1 


Cl 
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Regions 



4 a 

A 

B 

C 




k . 


4a 




M 






, ! 



Leading edge Trailing edge 




(a) Time interval I, 0«/c<l/(Af-K). 

(b) Time interval II, l/(M+l)<«/c<l/(M-l). 

(c) Time interval III, 1/(M— l)<t/c. 


Figure 4. — Variation of A p/g with chordwise station after sudden change in angle of 
attack. Regions defined in figure 3. 


that portion of the wing is unaware of the change in stream 
conditions. Over the forward portion of the airfoil (region 
A) the Ackeret-type loading corresponding to the modified 

c 

angle of attack is in evidence and for O’jjzrj extends over 

the entire chord of the wing. In region B the loading ex- 
periences a reduction in magnitude from the value over the 
forward portion of the wing. 






Leading edge Trailing edge 

(b) 


9 

Region 

4 to 0 ^ 

V 0 - 1 


and the loading has attained a static value agreeing with the 
Ackeret type of distribution. The loading in region B varies 
between the two constant values of regions A and C, dipping 
below that of region C and having its minimum value at 
,r=0. 

The loading produced upon entering a sharp-edged gust 
is pictured in figure 5. The loading in region C is zero, since 


Leading edge Trailing edge 

(c) 

(a) Time interval I, 0<t/c<l/(M+l). 

(b) Time interval II, l/(M+l)<t/c<l/(M-l). 

(c) Time interval III, l/(M—l)<t/c. 

Figure 5. — Variation of Ap/g with chordwise station after entering gust. Regions 
defined in figure 3. 
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The load distributions which have been developed were 
obtained for flight velocities in the supersonic regime. It is 
apparent, however, that the basic differential equation is not 
restricted to the case where M>1 and that the method of 
analysis affords a means whereby transient load effects can 
be studied for subsonic speeds. The essential difference 
between the latter problem and the results derived here lies 
in the relative position of leading- and trailing-edge traces 
in the x,t plane and the trace of the characteristic cone. 
Thus, for subsonic flight, the trace x= — t does not cut 
across the region occupied by the airfoil; whereas the cone 
stemming from the trailing-edge point c,0 does. A qualita- 
tive picture of the problem is obtained if the analogy between 
the nonsteady two-dimensional case and three-dimensional 
wing theory is used. The loading functions given in equa- 
tion (11) are equivalent to loading existing on a swept-for- 
ward tip of a three-dimensional wing. Thus, in figure 1 (a), 
x can represent distance measured spanwise, t can represent 
distance measured chordwise, and the shaded area can 
represent a portion of the plan form of the wing. Using this 
analogy, the loading which has been determined is merely 
load distribution over the swcpt-forward tip of a wing with 
constant chord and supersonic leading edge. When the 
case of the airfoil section traveling at subsonic speeds is to 
be considered, the problem becomes one of determining the 
loading over the swept-forward tip of a wing with constant 
chord and subsonic leading edge. 

DEVELOPMENT OF INDICIAL LIFT FUNCTIONS 

Since section lift coefficient c t is given by the expression 


c,=± (^dx 
CJ q 

the relations presented in equations (11) and (12) are suffi- 
cient for the determination of Ci at indicial lift coefficient for 
change in angle of attack, and c ts , indicial lift coefficient for 
an airfoil entering a gust. As a result of direct integration, 
the following results are obtained: 


First time interval 


0 <t< 


c 

I +M 



(13a) 


c 


4w 0 t 

cU„ 


(14a) 


c c 

Second time interval , ; - , j<t< 


4a 1 / 7r 

V- tt'L M\2 


+ arc sin 


1+M^^M-l 

c—Mt'' 


t 




arc cos 


Vm 2 - i 

4 w 0 t fir . . c—Mt\ 

2+ arcsin t ) 


t±Mc-tM 2 + l_ (c _ iM)2 " | (13b) 


J 


4 w 0 

arc cos 

ttFo VM 2 -1 


+ 

Afc+t 2 — MH 


(14b) 


Third time interval 


c 


<t 


4a 

l 


4w 0 

l ° UoVM 2 -! 


(13c) 

(14c) 


Values of the lift functions are plotted in figure 6 as func- 
tions of s, the distance traveled by the airfoil measured in 

half-chord lengths where «=“— The curves shown were 

calculated for values of M equal to 1.2, 1.31, and 1.46, since 



Figure 6.— Indicial lift functions ci a (s) and ci^(s) for various free-stream Mach numbers. 


the asymptotic values of Ci a and c t for the three cases agree 
with the values given in reference 1 for the subsonic wing at 
aspect ratios of °° , 6, and 3, respectively. No direct analogy, 
of course, can be made between the two cases. It is, how- 
ever, worthy of note that the variations in the indicial func- 
tions for the supersonic case are of the same order of magni- 
tude as those found in the finite-span incompressible case. 

From a knowledge of the lift function resulting from a 
sudden unit angle of attack, it is possible to express the lift 
corresponding to a given variable motion by considering 
the given motion as being composed of infinitesimal steps 
and summing the lifts corresponding to each step. Math- 
ematically, the problem corresponds to the use of the 
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Dtihamel or Faltung integral and leads to the so-called 
superposition theorem which can be written in the form 

c,(i')=«(0)c lo (i')+ J' Cjt'-r') dr' 

C la {t'^ T ') a (r')dr' (15) 

The primes on the variables in this equation indicate that 
true time is used. 


MOTION OF AIRFOIL IN GUST 

The results which have been obtained will now be applied 
to determine the forces on an unrestrained airfoil entering 
a gust. Since the motion of the wing is not prescribed it 
becomes necessary to equate the dynamical forces in order 
to relate the variables involved. Neglecting pitching 
moment and using Newton’s second law of motion, 

m forces (16) 

where w is the vertical velocity of the wing, m is the mass of 
the wing, and the forces to be summed result from the lift 
on the wing and the impressed force resulting from the 
action of the gust. By means of equation (15), equation 
(16) can be rewritten in the form 

m W + <lS df'fo C ^'-r')cc(r')dr'^ C Lg (t')q_S (17) 

Introducing a change of variables such that 

t i= T ~ a 0 

_ 2m 

the equation becomes, finally, 

M C*.(T~Tdc l (T l )dT l *=«± C Lg {T) (18) 

Since «(0) = 0, equation (18) can be integrated to give: 

1 r Cf ( r l)d r 1+ 1 r<7 (T-r.MTOdr^o ( 19 ) 

to Mjo 1 fj-j 0 

which is an integral equation of the second kind with a varia- 
ble upper limit. The solution to this equation can be obtained 
quite satisfactorily by means of Liouville’s method of suc- 
cessive substitution. Using the relation C L =n and per- 
forming the proper manipulations gives : 



and setting 

w 

V~ a > 


— 0 C L =C L (D-i [ T 0 Lu {T~l\)C Lg (7\)dT l + 

w 0 jujo “ £ 

±f*C L JT- TOdT t J* C La (Ti—T 2 ) C Lg (T 2 )dT 2 ~. 


( 20 ) 


Equation (20) is known to converge uniformly for 
and in the applications of this report the 


maximum lift was always experienced in the region of con- 
vergence. 

The values of lift coefficient C L determined from the solu- 
tion of equation (20) are shown in figures 7 (a), 7 (b), and 
7 (c) for various values of n and for M—1.2, 1.31, and 1.46. 
Figure 8 shows the variation of the maximum lift coefficient 
attained plotted as a function of the density parameter n for 
the same three values of M. Figure 9 furnishes a comparison 
between the maximum lift-coefficient increment given in 
reference 1 for an aspect ratio of 6 and the corresponding 



Figure 7.— Variation of increment of lift coefficient during flight through unit 
sharp-edged gust. 
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O 2 4 G 8 IO !2 14 

Half- chords traveled, s 

(b) Af=i.3i. 



(c) M=1.4G. 
Figure 7.— Concluded. 


value calculated in the present report for M=1.31. The 

4 772 

results are plotted as functions of *° correspond with 

the density parameter used by Jones. The correspondence 
which was noted for the indicial lift functions with Mach 
number replacing wing aspect ratio is still in evidence. 
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O 80 160 240 

Wing density parameter, fj.=2m/pScM 

Figure 8. — Maximum increment in lift coefficient attained during flight through unit 
sharp-edged gust for various values of free-stream Mach number. 



Figure 9.— Comparison of maximum increments in lift coefficients during flight through 
sharp-edged gust between subsonic wing arid supersonic wing. 

FORCES DEVELOPED ON GIVEN WING 

As an example of the uses to which the results just obtained 
can be applied, consider an airplane with wing of plan form 
such as the one shown in figure 2. Assume a wing loading 
of 40 pounds per square foot, a chord length of 8 feet, and 
let it be specified that the wing is flying at a Mach number 
equal to 1.2 and that the wing is built to withstand forces 
producing accelerations between —3 and 5 times gravita- 
tional acceleration. It is proposed to find at what altitudes 
the wing may be subjected safely to a gust possessing a 
vertical velocity of 50 feet per second. 

If F denotes the net lifting force on the wing, 

F=W+L g ( 21 ) 

where W is wing weight and L e is the total lift produced by 
the gust. If (A C L ) m ax is the maximum increment in lift 
coefficient attained in a unit gust and A is the acceleration 
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factor measured in multiples of g, then equation (21) is 
expressible in the form 

A=l + (AG L ) max ^-^ P0 V 0 
Since it is required that \A — 1|=4, it follows that 

{^C L ) max ^~ (22) 


Figure 10 shows the limit curve of (A C L ) max plotted as a 
function of flight altitude. From a knowledge of g—2 Wj 



Figure 10. — Safe altitudes for wing with loading of 40 Ib/ft 3 flying at M=1.20, based on 
50 ft/see gust velocity and maximum load factors of 5 g and —Zg. 


p 0 gMSc, however, the value of (A C L ) max actually attained 
by the wing entering' the gust can be calculated. Such 
values are also included in figure 10 and indicate that, under 
the given conditions, the wing should not fly at an altitude 
less than approximately 28,000 feet. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec. 5, 1947. 
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